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Abstract 



We study local and global approximations of smooth nets of curvature 
CN , lines and smooth conjugate nets by respective discrete nets (circular nets and 

planar quadrilateral nets) with edges of order e. It is shown that choosing 
the points of discrete nets on the smooth surface one can obtain the order 
(<~^ . approximation globally. Also a simple geometric construction for approximate 

I determination of principal directions of smooth surfaces is given. 

o 

> ■ 

: 1 Introduction 

^ Discrete conjugate nets, defined as mappings — > M'^ with the condition that 

each elementary quadrangle is flat, and discrete nets of curvature lines, defined as 
discrete conjugate nets with additional property of circularity of the four vertices of 
every elementary quadrangle, play an important role in the contemporary discrete 
differential geometry (see e.g. [2] for a review) and have important applications to 
computer-aided geometric design [T^ . 

This paper is devoted to a careful! study of the order of approximation of discrete 
nets to a given smooth conjugate net or the net of curvature lines on a smooth surface 
that can be achieved. We also present a remarkably simple approximate construction 
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of principal directions on a given smooth surface using the standard primitive of 
discrete differential geometry: an elementary circular quadrangle inscribed into the 
surface. Roughly speaking, our results show that the previously known upper bounds 
([3]) can be made one order better for conjugate nets and discrete nets of curvature 
lines. Moreover, one can impose one additional geometric requirement: all vertices 
of approximating discrete nets should lie on the original smooth surface. 

The paper is organized as follows. In Section [2] we present a method of construc- 
tion of principal directions in a neighborhood of a non-umbilic point on a smooth 
surface using only the four points of intersection of the surface with an infinitesimal 
circle. Our method is in fact based on a third-order approximation of the given 
smooth surface and gives an e^-approximation of principal directions if the infinites- 
imal circle has the radius e. Different methods for approximation of principal direc- 
tions have been proposed recently, see for example [S] . Re-meshing techniques based 
on determination of principal directions are widely used in contemporary computer 
graphics ([!]). 

Section [3] is devoted to the study of the order of approximation of a single elemen- 
tary quad of a discrete conjugate net (resp. discrete circular net) to the respective 
infinitesimal quadrangle of a the smooth net of conjugate lines (resp. curvature 
lines) on the original smooth surface. The results of this Section are used in Sec- 
tion m where we prove the main global results on the order of approximation of 
respective discrete nets to a given smooth conjugate net or the net of curvature 
lines on a smooth surface. 

We use the following terminology and notations: d{A, S) is the distance from a 
point A to another point 5* (or a line, plane, etc.). A point A is said to be e'^-close 
to S if d{A,S) = 0{e''), that is d{A,S) < Ce^ for e 0. A is said to lie at e^- 
distance from 5", if d{A^ S) ~ Ce^ for some C 7^ 0, we will denote this hereafter as 
d{A, S") ~ e^. The same terminology and notations will be used in characterization 
of the asymptotic behavior for other (numeric) quantities depending on e. 

2 Principal curvature directions from circles 

In this Section we show that one can use the elementary quadrangle of a discrete 
circular net inscribed into a smooth surface (see Section 13.21 for more details) for a 
rather precise construction of the principal directions on a smooth surface. 

• Euclidean construction. Take three points A, 5, C on a smooth surface / : 
— > without umbilic points, such that d{A, B) ~ e, d{B, C) ~ e and the 
angle ZABC is not e-close to or 27r.Then the circle u passing through A, B, 
C has radius of order 0(e) as well; from elementary topological considerations 
one concludes that u has at least one more point of intersection with /(fi). 
Denote it (or one of them) as D. Take the straight lines [AC), {BD) and the 
point of their intersection Z = (AC) fl (BD) (see Fig. [T]). The bisectors bi and 
62 of the angles formed at Z by (AC) and (BD) are obviously perpendicular. 
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In fact bi and 62 approximate the directions of smooth curvature hues on 
the surface: as we prove below, the directions of bi and 62 are e-close to the 
principal directions at any point M G f{^) which is e-close to A, B, C, D. In 
fact one can find the principal directions with only 0(e^)-error at some specific 
point Zf of f{^) which "lies over Z". Note that the order of the points A, B, 
C, D on the circle u is important for finding the principal direction with the 
error of order e^. 

Namely, let Zf be the intersection of /(fi) with the straight line / passing 
through Z and perpendicular to the plane Tr^ = (ABC). Take the orthogonal 
(parallel to I) projections pi and p2 of the bisectors bi and 62 onto the tangent 
plane tt/ of f{fi) at Zf. 

Theorem 1 For a general smooth surface /(fi) in a neighborhood of a non-umbilic 
point the Euclidean construction given above produces the directions pi, p2 which are 
e"^ -close to the exact principal directions at Zf. 

In order to prove this Theorem we establish at first some simple but remarkable 
facts about planar quadrics and Dupin cyclides (Lemma [2] and Theorem [3]) and give 
another (Mobius-invariant) construction of the principal directions (Fig. [2]). In fact, 
the intersection curve 6 = f{^) H tt^ is e^-close to a quadric — the Dupin indicatrix 
of f{^) at P (see the Appendix); this explains in particular why we may assume 
that the number of intersection points of 00 and f{^) is exactly four. If we take then 
an arbitrary circle a; in a plane, intersecting a quadric in points A, B, C, D (cf. 
Fig. [1] for the elliptic case; the order of the points A, B, C, D on the circle is not 
necessarily consecutive for Lemma[2]below to hold), then the directions of the angles 
formed by the straight lines {AC) and (BD) are parallel to axes of the quadric: 



Figure 1: Construction of the axes directions of a quadric: A, B, C, D are the four 
points of intersection with a circle, 61 and 62 are the bisectors of the angles between 
{AC) and {BD). The bisectors bi and 62 are parallel to the axes of the quadric 
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Lemma 2 For any plane quadric with four point of intersection with a circle the 
bisectors hi and 62 (shown on Fig. [I] /or the elliptic case) are parallel to the axes of 
the quadric. 

Proof. For simplicity we give the proof only for the case of a non-degenerate central 
quadric Q{x, y) = {Kix^ + — 2e^ = 0; the other cases can be proved along the 

same lines. The coordinates of the intersection points A, B, C, D satisfy the system 

r Q{x, y) = {Kix' + K^y^) - 2e' = 0, 

\ uj{x, y) = {x- x^y + {y - yef - = 0. ^ ' 

They also lie on the one-parametric sheaf of quadrics Q\: {Q{x, y) + \uj{x, y) = 0}. 
Obviously the directions of the axes of all Qx are the same. On the other hand, for 
some values of A this quadric is degenerate. These values are the three roots Aj of 
the equation 

Ki + X -2Axe 
K2 + X -2Xy, 
-2Xx, -2Xy, 2e^ - Xr] 

Two of them correspond to the degenerate quadrics composed of the pairs of straight 
lines {{AB),{CD)} and {{AD), {EC)] shown on Fig. [H The third degenerate 
quadric is the pair of straight lines {{AC)., {BD)}, shown on Fig. [TJ 

Since the bisectors of these pairs of lines define exactly the axes of the respective 
degenerate quadric Qx, we have proved that their directions are the same as for 
Qo = {Q{x,y) = 0}. □ 

Note that we have not used the ordering of the points A, B, C, D on the circle 
u or on the quadric, so the statement of Lemma [2] holds true also for the bisectors 
of the angles formed by the lines {AB) and {CD), as well as {AD) and {BC). On 
the other hand in order to obtain the approximation of order below we fix the 
ordering as shown on Fig. [H 

Next we need to reformulate the Euclidean construction given above in terms of 
Mobius geometry: 

• Mobius-invariant construction . We start with a smooth surface / : 
Q — > M^. Take a sphere S tangent to f{^) at some point P, such that 
for its radius R and the principal curvatures Ki, K2 {Ki < K2) of f{^) at 
P the inequalities Ki < 1/R < K2 hold (with non-infinitesimal differences 
\Ki — 1/R\ and \1/R — K2\ )■ In this case S intersects an e-patch on f{fl) 
along two smooth lines 71, 72 (Fig. [2]): Qf H S = jiU ^2, the angle between 71 
and 72 is non-zero and non-infinitesimal (this follows from the classical lemma 
of Morse, cf. for example [ID]). Take any circle u; C S* of radius ~ e with P 
inside such that the distances between P and the four points of intersection of 
Lj and Q^: ci;n7i = AUC, u;n72 = BUD are of order 0(e) as well. Construct 
two auxiliary circles Ui, UJ2 on 5* defined by the triples of points {A, P,C} and 
{B,P,D} respectively (not shown on Fig. [21 uji are very close to 7^). Then 



(2) 
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the bisectors 61, 62 of the angles formed by tui, uj^ at P obviously lie in the 
tangent plane to /(f^) at P. As we show below, hi and 62 are e^-close to the 
principal directions of /(f^) at P. 




Figure 2: Approximate Mobius-invariant construction of the principal directions. 
A sphere S is tangent to /(f^) at P, it intersects /(f^) along two smooth lines 71, 
72. A circle u; C 5* of radius ~ e with P inside gives four points: fl 71 = A U C, 
u;n72 = BUD, the points A, B, C, D should lie on e-distance from P. The bisectors 
bi, 62 of the angles formed by the circles Ui, uj2 (not shown) passing through the 
triples of points {A, P, C} and {P,P, D} are e^-close to the principal directions of 
f{9) at P 

First we prove that for Dupin cyclides (standard Mobius primitives, playing the 
role of the osculating paraboloids of the classical Euclidean differential geometry, 
see OE]) &i are exactly the principal directions: 

Theorem 3 For a smooth non-umhilic point P on a Dupin cyclide D, the bisectors 
bi and 62 of the Mobius-invariant construction give the principal directions of T) at 
P. 

Proof. Since our construction of Ui is Mobius-invariant, we may transform D into 
one of the normalized Dupin cyclides: a torus, a circular cone or a circular cylinder. 
We give below a detailed proof for the case of a torus; the other cases can be easily 
proved along the same lines. 

Making, if necessary, another Mobius transformation, we can reduce the config- 
uration to the following: the torus is given by the equation 

the point P has the coordinates (P+r, 0, 0) and the sphere S is given by the equation 

{x + p- {R + r)Y + + z'^ = p'^. (4) 
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After a Mobius inversion [x, y, z) — > (xi, yi, ^i) with the center P: 
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(5) 



the sphere ([1]) will be transformed into a plane xi = const = fi. We prove now that 
the intersection of the image of the torus with any plane xi = is a quadric of the 
form ayf + bzf = k. In fact, from ([5]) we have: 



so 



X - (i? + r)Y + y^ + z^ 
[X — K — r] / u — I (x 



R-r)/fi- 
[I3l) may be now reshaped to 

(^(x -{R + r)) /fi + 2{R + r)x 



= [x — R — r)/ ji 



(6) 
(7) 



2rR 



= AR\x^ + y^), 

so one can express y"^ in terms of x, R, r, fi. Substituting ([7]) and the found 

expression for y"^ into ayf + bzf = {ay"^ + bz'^) j {^{x — {R + r)Y + y"^ + z^^ one 

gets a rational expression containing only x, the constants i?, r, /i and a, b. A 
straightforward calculation shows that for a = r{l + 2fi{R + r)), b = (1 + 2r/i)(r + 
this expression is a constant n = —R{1 + 2/ii?)(l + 2/ir)(l + 2/ir + 2/i_R)/4. 

Since the images of the circles cui, u!2 defined above are now straight lines pass- 
ing through the points of intersection of the image of the circle u; and the found 
quadric on the plane Xi = fi, the statement of this Theorem is now equivalent to 
the statement of Lemma [2J □ 



Theorem 4 For an arbitrary smooth surface f{^) in a neighborhood of a non- 
umbilic point P the Mobius-invariant construction produces the directions bi, 62 
which are e^-close to the exact principal directions at P. 

Proof. First we perform a Mobius transformation that maps the sphere 5* and the 
circle u into themselves and brings P into the (Euclidean) center of u on S, i.e. into 
the point P' e 5* such that the distance from P' to all the points of u are equal. 
One can easily check that our requirement d{P, A) ~ e, . . . , d{P, D) ^ e guarantees 
that after the transformation these distances will remain of order 2^ e. Inside the 
circle a; on S* and in its e-neighborhood in where f{^le) lies, the Jacobian of this 
transformation is limited so the distances of order O(e^), p >2 will be transformed 
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into distances of the same order. Since the angles do not change, e^-close directions 
at P remain e^-close at P' and vice versa. Introduce now a Cartesian coordinate 
system such that P' is its origin, the coordinate plane {xy) is tangent to S (so to 
/(fi) as well) and the x, y-ax.es are the principal directions of f{^) at P' . The plane 
where uj lies will be given by the equation z = Cie^, the circle u will be defined 
as + = c^e^ and fX^) will be defined in a neighborhood of P' by its Taylor 
expansion 

z = {Kix^ + K2y^)/2 + (agox^ + a2ix% + oizxy^ + ao^y^) + o(e^). (8) 

So the intersection points Pi = A, P2 = B, P^ = C , P^ = D oi uj and f{^) will 
have respectively the coordinates Xi = exio + e^xn + o(e^) with (xjo, yio) being the 
solutions of 

X + y = C2, 
(i^ix2 + AV) = 2ci, 

so we see that 

(a;io, 2/10) = (2:20, -y2o) = (-a^so, -yso) = (-3^40, 2/4o)- (10) 
The next terms {xn^yn) are found after substitution of {xi,yi) into 



x^ + y^ = 4, 



'11] 



{Kix^ + K2y^)/2 + e(a3ox3 + aaix^y + ai2xy^ + aos^^) = Ci, 
cancellation of the second-order terms using (Q and retaining only third-order terms: 
2xioXii + 2yioyii = 0, 

KiXioXii + K2yioyii + (osoxfo + a2ia;?o2/io + ai2Xiol/I) + (^osvfo) = 0- 
From (fT2l) and (fTOj) we easily conclude that 

(2^11, 1/11) = (3:31, Z/31), (a^2i, 2/21) = (a;4i, Z/4i), (13) 

so the 0(e^)-shifted lines In, I21 defined by the pairs of points {(exio + e^xn, eyio + 
e^yii), {ex3o + e^Xsi,ey3o + e^y3i)} and {(exzo + e^^ai, e?/2o + e^y2i), (6X40 + 6^X41, e?/4o + 
6^1/41)} are parallel to the lines /lo, /20 defined by the pairs {(6X10, ^yio), (6X30, ^yso)} 
and {(6X20, 6^/20), (6X40, 6|/4o)}- 

Thus the bisector directions between In, /21 in the plane are the same as for the 
unperturbed lines Zio, /20; so they are obviously parallel to the coordinate axes — 
the principal directions of f{^) at P' which is e^-close to both Zq = Iiq fl /20 and 
Zi = III n /21. 

Now we have two pairs of circles in M^: uio, U20 passing through the triples 
of points {(6X10, eyio),P', (exso, ey^o)} and {(6X20, mo), P' , (6X40, 6I/40)} and the pair 
ujii, UJ21 passing through {(6X10 + 6^X11, eyio + e^yn), P' , (exso + 6^X31, e?/3o + e^y-n)] 
and { (6X20 + 6^X21 , 6?/2o + 6^?/2i ) 5 -f"; (6X40 + 6^X41 , e?/4o + 6^?/4i ) } . One easily concludes 
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from (fT3l) that ujiq, uju are tangent at P', as well as UJ20, ^21- Thus their bisectors 
are the same and parallel to the axes; the next 0(e^)-terms in approximations for 
the intersection points A, B, C, D will give a 0(e^)-change in the directions of bi, 
b2. □ 

Remark. From the proof of the previous Theorem we see that in fact the Mobius- 
invariant construction gives the same (with O(e^)-error) principal directions if we will 
substitute instead of the original surface f{^) its tangent paraboloid at P (after the 
Mobius transformation described in the beginning of the proof). Also it should be 
noted that the 0(e^)-approximation was achieved taking into account approximating 
paraboloid of third order; the results obtained in this Section are therefore third order 
results despite the seemingly second-order construction based on an infinitesimal 
circle. 

Proof of Theorem [H Now we are in a position to prove the statement about 
0(e^)-approximation of principal directions at Zj in our Euclidean construction. 
For this we first perform an auxiliary construction in M.^ (Fig. [3]): 



C 




D 



Figure 3: Auxiliary construction 

take the circles uji, UJ2 defined by the point triples {A, Zj, C} and {B, Zf, D}, then 
for the sphere S where all these points and circles lie we take its tangent plane tts 
at Zj. Next we construct another two circles (i^ as the sections of the sphere 
S by the "bisectoral planes" passing through Zf and the lines 61, 62 respectively. 
First we remark that the constructed plane vts is e^-close to the tangent plane vry 
of /(f2) at Zf. the circles cui, 1^2 are approximating the osculating circles of the 
planar sections of ji^l) by the planes (A, Zj, C) and (-B, Zj, D\ As one can check 
by a direct computation for any planar smooth curve 7, if one takes the circle lo 
passing through three points lying on e-distances from each other on 7, then the 
tangents to lo and 7 at any of the three pins will be e^-close. So the tangents to 
the aforementioned planar sections and the tangents to the respective circles cji, 102 
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at Z/ are e^-close. (Here we use the same terminology of "e^-closeness" for pairs of 
lines or planes, this means that the respective angles are 0(e^).) This shows that 
our sphere S is "very close" to the sphere S (which should be tangent to f{^) at 
Zf = P) used in the Mobius-invariant construction. We also see that our Euclidean 
construction is nearly reduced to the Mobius-invariant construction, with the sphere 
5" substituted by S: the orthogonal projections (along the line / = (ZZf) ± txabcd) 
of hi, 62 onto TTs are nothing but the tangents to the circles (32 which are e^-close 
to the circles which exactly bisect the angles between the circles a;i, a'2 at Zf. This 
follows from the fact, that the angle between tiabcd and tt/ is 0(e): as one can 
easily check by a direct computation, in this situation the (non-infinitesimal) angle 
between any two lines on one of the planes and the angle between their orthogonal 
projections on the other plane are e^-close. 

Let us now estimate the angles between the tangents to the circles (lying in 
TTs) and the principal directions of /(fi) at Zj (lying in tt/, which is e^-close to tts). 
For this we use the same technology as in the proof of Theorem HI first we perform 
a Mobius transformation leaving S and uo invariant and bringing the point Zf to 
the center O of on E, then introduce a Cartesian coordinate system with O as its 
origin and vr^ being its (x?/)-plane. Now the equation ([8]) of the surface f{^) will 
be modified by small linear terms: 

z = {Kix^ + K2y'^)/2 + (asox^ + a2ix^y + ai2xy'^ + aosZ/^) + crix + (J2y + o(e^), (14) 

with (Tj = O(e^). The first terms Xjo, Vio in the Taylor expansions of the coordinates 
of the points A, B, C, D will be found from the same system iQ so we have the 
relations ffTOj) . The next terms xn, yn are found from a modified system of the form 

m-- 

{ 2xioXii + 2yioyii = 0, 

\ KiXioXii + K2yioyii + (asoa^fo + «2ia;fol/io + auXtoy^Q + ao^yf^) + ^x-io + ^yio = 0. 

We see that the main conclusion (|T3i) still holds; from this point we just follow 
the guidelines of the proof of Theorem H] and show that the bisectors of the angles 
between the circles ui, 002 coincide (up to a O(e^)-error) with the directions of 
the X- and y-axes. On the other hand since we know that cTj = O(e^), using the 
standard differential-geometric formulas for calculation of the coefficients of the 
second fundamental form and the principal directions we conclude that the principal 
directions of f{fl) given by the equations ([H]) are e^-close to the x- and y-a.xes at 
Zf as well (as the directions of the straight lines in M.^). □ 

Remark. From the proof of this Theorem we see that in fact one can assume 
(with O(e^)-error) in the Euclidean construction that the tangent plane at Zf is the 
readily constructible plane tts. 
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3 Local results 



3.1 Smooth and discrete conjugate quads 

Suppose that a smooth surface f{^) parametrized locally by some curvilinear net of 
conjugate lines is given: / : Q — > M.^, C = {{u,v)}. Take some initial point 
A = f{uo,Vo) and points B = f{uo + e,Vo), C = /{uq^vq + e) at e-distance on the 
two conjugate lines of the net on f{fl) and let D = f{uQ + e,fo + e) be the fourth 
point of the curvilinear quad on f{VL). 
Using the conjugacy condition 

fuv = a{u,v)fu + b{u,v)fy (15) 

and its derivatives, one can easily estimate the distance from this fourth point to 
the plane ttabc defined hj A, B, C; we formulate this as our next Theorem. 

We remind that smooth conjugate nets are supposed to be non- degenerate, that 
is the directions of the curvilinear coordinate lines of such a net are non-asymptotic, 
so the angle between the tangent vectors and fy is everywhere different from zero. 

Theorem 5 For an arbitrary smooth non-degenerate conjugate net f{^), d{D, ttabc) 

0(6^). 

Proof. Using the standard Taylor expansions one gets: 

AB = efu + y/uM + ^fuuu + o{e^), 

AC = efy + f + + o(e3), (16) 

AD = AB + AC + e^fu. + ^(/™„ + fuuv) + o{e^), 

where all derivatives are taken at the point {uq,Vq). 
Taking into account (fT5|) and its derivatives 

fuuv (^ufu ~l~ byfy -\- Ojfy^u -\- bijifu ~\~ bfy^, 

fuuv (^vfu ~l~ byfy -\- (X(Clfu -\- bfy) -\- bfyy, 

we can compute the triple product {{AB x AC) ■ AD) which gives the volume of the 
parallelepiped spanned by the vectors AB, AC, AD: 

{{AB xAC)-AD) = e\AB x AC) ■ fyy + ^(AB x AC) ■ (/„™ + fuuv) + o(e^) = 



2 



2 

{fu ^ fv) ' \ fuvv ~\~ fuuv ) ~l~ {fuu ^ fv) ' fuv ~\~ {fu ^ fvv) ' fu 



+ 0(6^1 



e 



a{fu X fy) ■ fuu + b{fu X fy) ■ fyy + a{fuu X fy) " fy + b{fu X /„„ ) " /„ + o(e^) = o(e^) 
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Since the area of the base of this parallelepiped \AB x AC\ ~ e^, we get that its 
height d = o(e^), so at least we may state that d = 0{e*). In fact, using Taylor 
expansions in (fT6|) up to order 0{e'^) one can obtain after a lengthy computation 
that 



( {AB xAC)-AD) = — [(a, - ab) ( x /„)■/, - (6, - ab) x /J ■ /„J + o(e^) (17) 

which proves that for a generic conjugate net we have d ~ e^. □ 

In fact we can even choose M (the fourth point of an elementary planar quad 
ABCM) on the plane ttabc sufficiently close to D and lying on the given smooth 
surface f{fl): 

Theorem 6 For a given smooth non- degenerate conjugate net f : Q — > M and the 
plane ttabc constructed as above one can choose a point M G habc H /(fi) such that 
d{M,D) = 0{e^). 

Proof. As we prove in the Appendix, the intersection of the plane ttabc and the 
surface f{^) in the e-neighborhood of the points A, B, C, D for sufficiently small 
e is a curve Iabc close to a quadric - the Dupin indicatrix of fX^), and the angle 
between tiabc and /(fi) is ^ e in all points of Iabc- According to Theorem [5], 
d{D,7rABc) = 0{e^). From this we can easily see that the distance between D and 
the closest to D point M on Iabc should be 0(e^/e) = O(e^). □ 



3.2 Curvilinear quads of curvature lines and plane circular 
quads 

For a given smooth surface f{^) parametrized by curvature lines one can prove a 
similar result, if we take the circle ujabc passing through the points A, B, C defined 
as in Section 13.11 

Theorem 7 For arbitrary smooth net of curvature lines in a neighborhood of a non- 
umbilic point A on f{fl), the distance between the fourth point M of intersection of 
the circle ujabc with f{fl) and the point D has order 3: d{D, M) = O(e^). 

In order to prove this theorem we need to establish a few auxiliary Lemmas. 

Lemma 8 For arbitrary smooth net of curvature lines in a neighborhood of a non- 
umbilic point A on fiVL), d{D,uJABc) = 0{e^). 

Proof. Following an approach proposed in jlj, we introduce for each of the points 
A, B, C, D purely imaginary quaternion A, B, C, D. The point A will be chosen 
as the origin, i.e. A = 0. The basic quaternions I, J, K are chosen to be tangent 
to the curvature line directions /^t, f^ and the normal vector to the surface in the 
initial point A respectively. In view of future applications in the theory of triply 
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orthogonal coordinate systems we include the given surface f{^) into such a system 
f{u, V, w), f{fl) being one of the coordinate surfaces f{u, v, w = Wq). This is always 
possible (see e.g. [7]): one can take for example the one-parametric family of surfaces 
parallel to f{^) and the other two one-parametric families of developable surfaces 
defined by the curvature lines of f{^) and the normals to /(fi). 

Introducing for this 3-orthogonal system the Lame coefficients Hi{u,v,w) = 
\dif\, di = d/du, 82 = d/dv, = d/dw, normalized vectors Vi = dif /Hi and the 
rotation coefficients (3ik{u) = diH^/Hi, i ^ k, (3ii{u) = 0, we have the following 
relations ([7j): 

diHk = l3ikHi, 
diVk = PkiVi, 

<9^^^ = -Es^./3s^K, (18) 
djPik = PijPjk, j ^ k, 

where i,j,k G {1,2,3}, i ^ j ^ k. In the initial point A we have Vi{uq,Vq) = I, 
V2(uo, vq) = J, Vs^uq, vq) = K. The Taylor expansions for the other points are now 
easily obtained after differentiation of /„ = HiVi, = H2V2 using (ITSl) : 

B = eH^l + ^{d.H^l - ^2iH^3 - f3siH,K) + o{e''), 

C = €H2l + f (-/5i2i/2J + d2H2J - /332^2K) + o{e^), (19) 

D = B + C + e^{p2iH2l - Pi2Hi3) + o{e^). 

Calculating now the quaternionic cross-ratio (|1]) Q = (A— B)(B— C)~^(C— D)(D — 
A)~^ and taking its imaginary part, one can see that ImQ = o(e). According to [1] 
we conclude that p = o(e^) for the distance p from the point D to the circle ujabc 
of size ~ e defined by the points A, B, C. □ 

In the following we fix an e-patch /(fi^) of size ~ e on the surface f{^), such 
that all the points A, B, C, D belong to it. All subsequent constructions will be 
applied only to this e-patch /(l^e)- 

From Lemma [16] (proved in the Appendix) we know that the intersection curve 
I ABC = /(^) nTTABc IS e^-close to the Dupin indicatrix of f{VL) at some point M. In 
fact we can take the indicatrix at the following point P instead: choose the "center 
point" P on /(f2) such that P is the intersection of f{^) and the normal to the plane 
T^ABC passing through the center of the circle ujabc- In the Cartesian coordinate 
system (similar to the system used in the proof of Lemma [T6|) . where the osculating 
paraboloid at P is T = {2; = {Kix^ + K2y'^)/2}, the points A, B, C, D will have 
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coordinates 



A(-^ + o(e), -f + o(e), M2e2 + o{e^: 

BU + o{e),-^,+o{e),^e' + o{e': 
c(-f + o(e),f + o(e),Mae2 + o(e2; 

D(f + o(e),f + o(e),M.e2 + o(e2; 



(20) 



so for the angle 6 between the plane ttabc and the tangent plane vrp to /(f^e) at P 
we have 6' = O(e^). Moreover, the following is true: 

Lemma 9 The centers of the Dupin indicatrix Qp = 7 {~\ habc o-nd ujabc ore en- 
close. The angles of intersection of uoabc o-nd Iabc (if^^ ~ e'^ = 1. 



Proof. From ( 1201) we easily deduce the first statement. The radius of uabc is 
eV2 + o(e), while the axes of Qp are ^^^^ + o(e) and ^^^^ + o(e). Since we 
assume \K2 — Ki\ y e on Q we see that the intersection angles of Qp and uabc, so 
also of uabc and Iabc, are ~ e° = 1. □ 

Now using the same technique as in the proof of Lemma [121 we deduce that the 
Dupin indicatrix at P and Iabc are e^-close. 

Proof of Theorem\^ We know (Theorem [5]) that D is e^-close to ttabc- Using 
the fact that the angle of intersection of habc and /(fij is ~ e (so >- e^) we see that 
D is e^-close to Iabc = t^abc H f{fl)- On the other hand (Lemma [8]) D is e^-close 
to the circle ujabc- Since the angles of intersection of ujabc and Iabc are ^ e° = 1, 
we see that D and M = Iabc H ^abc are e^-close as well. □ 



4 Global results 
4.1 Conjugate nets 

Using the results of Section [3TT] one can try to construct inductively an approximating 
discrete conjugate net for sufficiently small e > 0. At the first glance the following 
simplest strategy may be applied: starting from the initial point A = f^Q = /(mq, Vq) 
on a finite piece of a smooth surface / : Q — > M.^ parametrized by conjugate 
coordinate lines, fix two series of points fi^ = f{uo + ie,vo), /o,i = f{uo,vo + ie) at 
e-distances on two curvilinear coordinate lines on f{VL) passing through A. We can 
construct the next point fl^ as the orthogonal projection of /n = /(mq+c, VQ+e) onto 
the plane passing through /qo, /lo, /oi! then f^i as the orthogonal projection of /21 = 
/(Mo + 2e, fo + e) onto the plane passing through /lo, /20 and fl^, etc. This inductive 
process is shown on Figure H] with initial points fij = f{uQ + ie,vo + je) on the 
surface f{^) with smooth coordinate lines (shown as dash-lines). The approximating 
discrete conjugate net consists of the points //j, //g = fio, foi = foi- 

This "geometric analogue" of the analytic approximation results of [31 [H] should 
give us, if the results of [3l[TT] were directly applicable, an estimate c?(//,-, fij) < C{i+ 
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/so = /so 



Figure 4: A discrete conjugate net f^j and a smooth conjugate net fij which are 
e^-close. The points of f[j and fij on the initial curves coincide 

j)-e^, so for the complete simple piece f{fi) we would have d{flj, fij) = O(e^). In fact 
d{fij, fij) shows quadratic behavior: d{f[j, fij) < Ci-j-e^, so on fip) we have a much 
weaker estimate d{flp fij) = O(e^). Indeed, a straightforward calculation shows that 
for the distances fij) of the four vertices of any elementary discrete quadrangle 
one has d{f^'^^j^^, /i+ij+i) - f^+l,j) + d{flj+i, f^,j+l) - d{flj, f,j) + e^(j){uQ + 

ie, Vo + je) where v) = (a„ - a6)(/„„ x /„) ■ - (6„ - a6)(/^,^, x /„) ■ /„ is the 
coefficient of in the right hand side of f|T7j) . Then for example for special surfaces 
which have constant (f){u, v) = (po one easily obtains d{flj, fij) ~ z ■ j' • e"*. 

Paradoxically, a better strategy consists in choosing the points f-j of the approxi- 
mating discrete net on the given smooth surface /, although Theorem [6] suggests that 
one has worse local error ^ than the local error ~ of Theorem [51 The construc- 
tion we propose below is applicable to arbitrary simple pieces f{^) without parabolic 
points, i.e. points where along with the diagonal coefficient {f^v x fu) ■ fv of the sec- 
ond fundamental form (cf. (fTSl) ) at least one of the other coefficients x /„) ■ 
ifvv X /„) ■ f\, also vanishes. 

First we fix the vertices //g = fio = f{uQ + ie, ^o) and all even "columns" /|j j = 
f2i,j = f{uo + 2ie, Vo + je) in the new approximating discrete conjugate net (marked 
on Figure [5] as (^-points) . 

Then we proceed as follows. Two planes defined by the triples of points {foo, fio, foi}, 
{/lo, /20) /21} intersect along a straight line passing through /(q = fiQ. As we show 
below in Lemma [TU|, this straight line intersects f{^) in another point f^^ which 
is close to /ii, = 0{e^), whereas the distance between this new point 

fl^ and the line v = Vq + e of the smooth conjugate net is of order O(e^). Other 
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Figure 5: A discrete conjugate net //^ with its vertices on f{^), which are en- 
close to the vertices of a smooth conjugate net fij. The original smooth net fij is 
schematically shown as a rectangular plane net 

points f^i, f^i, . . . of the first row are constructed in the same way and obviously 
have the same error estimates. The second row /2j+i 25 = 0,1,... of the vertices 
of the approximating discrete net is constructed analogously using the fixed points 
/2fc2' ^ = 0, 1, ... of the second row on f{fl) and the constructed complete set of 
the points /,f^ of the first row (cf. Figure [5]). 

As Lemma [TT] below shows, /^j+i 2 ^^e now e^-close to the smooth line v = vo + 2e 
whereas the distance between 2 and f 21+1,2 along this line (the difference of their 
u-coordinates) is doubled: 5u{fl2) = ^u{fl]) + ii{uo,vo + e)e^ + o{e^) = 2fi{uQ,vo)e^ + 
o{e'), Svifl,) = -Svin,) + u{uo,vo + e)e' + o{e^) = o{e'), 

Here and afterwards we use the following notations: Sv{f^j) = f (//j) — v{fij), 
^"^ifij) = ^(/ij) ~ "^{fij)^ where u{P) and v{P) denote the respective curvilinear 
coordinates u, v for any point P on the surface /(fi). 

For the third row we will have 5v{fl.^) = 3/i(Mo, fo)e^+o(e^), 6u{fl.^) = i^(mo, vo)e'^+ 

This linear "w-drift" of the points flf. in contrast to " w-oscillation" of the same 
points is easily explained: the position of the points /fo = /io = /(^o + i^iVo) 
on the initial conjugate line v = vq may be changed using a reparametrization 
u I — > u = f{u) of the first conjugate coordinate u; thus the u-coordinates of the 
points /j\ of the first row may be considered as some new choice of u-parametrization; 
all subsequent rows simply shift accordingly. 

Lemma 10 For two adjacent infinitesimal curvilinear coordinate quadrangles ABDC 
and BDFE of e- size (see Fig. on a smooth conjugate net f : Q — > M.^ without 
parabolic points one can find in a -neighborhood of D a unique point M on f{^) 
such that the quadrangles ABMC and BMFE are planar. The curvilinear coordi- 
nates of M are u{M) = uq + e + Sum, v{M) = vo + e + Svm, with 6um ~ /^(%, ^^o)^'^; 
6vm ~ u{uo,vo)e'^. 
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Proof. We will follow the guidelines of the proof of Theorem [51 Using the Taylor 
expansions up to order O(e^) for the Cartesian coordinates of A{uq, vq), B{uo + e, Vq), 
C{uq,Vo + e), D{uo + e,VQ + e), E{uo + 2e,fo), F{uo + 2e,V(, + e) and some point 
Di{uo + e + Sui,Vo + e + 6vi) with some 6ui = O(e^), 6vi = O(e^) and taking 
into consideration fll5l) . we compute the triple products Wi = {{AB x AC) ■ ADi) = 

e^(no,t;o) + f(^Wi(/n.x/J-/,+5t;i(/,,x/J-/,)+o(e6),1^2 = {{BE xBd[) -FB) = 

e^'y{uo,VQ) + x Q ■ U - Svi{fy^ x /„) ■ /^j + o(e^) where 'j{u,v) is a 

combination of the coefficients f{u, v), l3{u, v) of ( fT5l) . their derivatives and the triple 
products (/tin X ■ /„, (/„, X /„) ■ fy. So if both of the latter triple products do 
not vanish, one can choose unique 6ui = /i(uo, vo)e^, 6vi = such that Wi = 0(e''), 

W2 = 0(6^). 

This shows (as in the proof of Theorem [6]) that Di lies in e^- neighborhood of the 
curves Iabc: Ibef of intersection of /(fi) with the planes (ABC), (BEF). Since 
the angle of intersection of Iabc: -^bef is — 1, we conclude that the point M of their 
intersection, close to Di, is in fact e^-close to Di. 

Using Taylor expansions for the same points up to order O(e^), one can prove 
that for a generic surface actually 6um ~ fJ'{uo,vo)e^, 5vm ~ z/(Mo,fo)e^ with /i ^ 0, 
z/ ^ 0. □ 




A = foo B = fio E = /20 



Figure 6: The adjacent infinitesimal smooth quads are defined by the points A, B, 
C, D, E and F (schematically shown as rectangles) and the planar quads ABMC, 
BMFE 

If one repeats the same calculation using instead of 5 a point Bi(uq + e + 
6u{Bi),vo + 6v{Bi)) with some 6u{Bi) = O(e^), 6v{Bi) = O(e^), one will get from 
the conditions Wi = 0{e'^), W2 = O(e^) that Sui = 5u{Bi) + /i(uo,fo)e^ + o(e^), 
5vi = -5v{Bi) + iy{uo,vo)e'^ + X{uo,vo)6v{Bi)e + 'y{uo,vo)5u{Bi)e^ + o(e^) where 
/i(no, Vq), ^'{uq, Vq), X^Uq, Vq) , 7(^0, fo) are some algebraic combinations of the coeffi- 
cients /(m, v), P{u, v) of (fT5|) . their derivatives and the triple products (/„„ x fu) ■ fv, 
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ifvv X fu) ■ fv Hence the following estimates hold for the second point M of inter- 
section of IabiC and IbiEF- 

Lemma 11 For two adjacent infinitesimal curvilinear coordinate quadrangles ABDC 
and BDFE ofe-size on a smooth conjugate net f : Q — > M.^ without parabolic points 
and some point Bi{uq + e + 5u{Bi), vq + 5v{Bi)), 5u{Bi) = 0{e^), 6v{Bi) = 0{e^), 
one can find in a e"^ -neighborhood of D a unique point M on /(fi) such that the 
quadrangles ABiMC and BiMFE are planar. The curvilinear coordinates of M are 
u{M) = UQ + e + 5uM, v{M) = VQ + e + 5vM, with 6um = Su{Bi) + fi{uo, vo)e^ + o{e^) , 
Svm = -Sv{Bi) + z/(uo,fo)e^ + X{uo,vo)6v{Bi)e + -f{uo,vo)6u{Bi)e'^ + o{e'^). 

Remark. As before, one can explain the "e^-tolerance" along the w-lines simply 
by removing this w-shift using a reparametrization of the curvilinear coordinate u. 
Note that we need this reparametrization only locally, for a given pair of elemen- 
tary infinitesimal quadrangles, in order to establish this estimate for them, and not 
globally, for the complete discrete conjugate net. 

For global estimates we need the following simple result ([3]): 

Lemma 12 (Discrete Gronwall estimate) Assume that a nonnegative function A : 
N — > M satisfies 

A{n + 1) <{l + eK)A{n) + K (21) 
with nonnegative constants K , k. Then 

A{n) < (A(0) + nn) exp{Kne). (22) 

Theorem 13 For a smooth conjugate net without parabolic points f : Q ^ M.^ and 
sufficiently small e > 0, there exists a discrete conjugate net with all its points 
on f{n), such that d{f,„ /^.) = 0{e^). 

Proof. From Lemma [11] we immediately see that for each column k = 2i + l the func- 
tion A{n) = \Su{f^^) \ satisfies the estimate fl2Tl) with K = 0, k = 2e^ ma.X(^u,v)en |a*(w, v)\ 
for sufficiently small e, so for all n we have (!22!) . Since the number of steps in each 
column is ^ 1/e we obtain the global estimate |5M(/2i+i.j)| < Ce^- 

Now we set A(n) = \^v{f^^)\ for each column k = 2i + 1. Taking the estimates 
of Lemma [TT] and the already obtained global estimate for |5M(/2j+ij)| we get (l2T|l 
with K = max(^u,v)en |A(M,f)|, k = 2e'^max(„^^)gn(|/i(M, w)| + \j{u,v)\), so ([22]) gives 
the global estimate |5^'(/|^+lJ•)| < Ce^. □ 

Remark. As we have observed above, the w-shift of the points f2i+ij has a 
much more interesting oscillating behavior; in particular one may conjecture that 
= 0{e') globally. 
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4.2 Curvature lines and circular nets 



We start the construction of a circular discrete approximation of a given finite piece 
/(n) of a smooth surface without umbihc points parametrized by curvature hues 
fixing points f-^ = fio = /(mq + ie,vo), f^^ = foi = f{uo,vo + ie) on two initial 
curvature hues passing throw a point /oo = f{uo,vo) (cf. Figure [7]). 




Figure 7: A circular net Z/^- with its vertices on /(fi), which are e^-close to the 
vertices of the smooth net of conjugate lines fij 

The first circle is defined by the triple {/oo, fio, foi}, its fourth point of intersec- 
tion with the surface will be chosen as the next constructed point fii of the discrete 
circular net. On the next step we can define two circles by the triples of points 
{/oi, /o2, /ii} and {/lo, /205 /n}, the fourth points of intersection of these circles 
with the surface will give us /li- The behavior of the new points on the 

third step is of crucial importance. Their error estimates can be obtained from the 
following generalization of Theorem [71 

Lemma 14 Let the points A, B and C on f{fl) be chosen, such that u{A) = uq + 
6ua, v{A) = vq + 6va, u{B) = uq + e + 6ub, v{B) = vq + 6vb, u{C) = uq + 6uc, 
v{C) = vo+e+Svc , with allSu, 6v of order = O(e^). Then the curvilinear coordinates 
of the fourth point M of intersection of the circle ujabc with f{Q) are: u{M) = 
u{fu) + Sum, v{M) = v{fu) + 6vm, with 

Sum = Sua + Sub - Suc + e{SvB - SvA)fii{uo, vq) 

+e{Suc - SuA)fJ.2{Uo, Vq) + /i3(%, Vo)e^ + yU4(%, Vo)e^ + o(e^), -^gx 

Svm = Sva - Svb + Svc + e{SvB - SvA)jyi{uo, vq) 

+e{Suc - Sua)i^2{uo, vq) + 1/3(^0, ^'o)e^ + 1^4(^0, vo)e'^ + o^e"^), 
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Proof. We will use the same quaternionic cross-ratio Q = (A — B)(B — C)~^(C — 
D)(D — A)~^. The obvious changes are necessary in the Taylor expansions flTIJl) . 
where we shall take into consideration the 6u- and 5u-shifts of the points; the ex- 
pansions themselves should include terms up to order O(e^). For simplicity we will 
assume 6ua = 0, 6va = 0, the general case requires only a shift by 6ua, Sva in the 
net of curvature lines on f{Q). 

Fixing for the moment the points A, B, C and introducing D G f{^), u{D) = 
uq + e + Sud, v{D) = vq + e + 6vd with some 6ud = 0{e^), dvo = 0{e^), one can 
compute the imaginary part of Q: 

/m(Q) = l(pi(Mo,^^o)e^ + 0H2{5vd - 5vc + 5vb) + o{e^)) 

-f j(p2(Mo,^^o)e^ + 9Hi{-5ud - Suc + Sub) + o{e^)) 

+K(p3e2 + e-^H^H2{Hl{5uB - 5ud - 5uc) + HI{5vd + Svc - Svb)) + o{e^)), 

where 9 = \{HlH2fiz2 - HiHlP^i) = HfHl{K2 - Ki)/2, Ki being the principal 
curvatures at the point A. Since we assume Ki ^ K2, choosing 

Sud = Sub - Suc + e^Hsiuo, Vq) + o(e^), 
Svo = -5vb + Svc + 6^1^3(^*0, vq) + o{e^), 

one will achieve Jm(Q) = o{e^)l + o(e^)J + o(e^)K. As a lengthier computation 
shows, adding appropriate fourth-order corrections in Taylor expansions, given in 
(|23|) (set there 6ua = 0, 6va = 0), we get /m(Q) = o(e^)I + o(e^)J + o{e^)K. This 
means that the point D with 6u£), 6vd given in fl23|) is e^-close to the plane {ABC) 
(this can be also checked using the triple product {AB x AC) ■ AD) and e^-close to 
the circle uabc passing through A, B, C. Using the the fact that the angle between 
/(fi) and the plane (ABC) is ~ e at D, we conclude that for the fourth point M of 
intersection of uabc and f{fl) we can keep the same expansions □ 

Remark. Linear behavior of 6ud and 6vd in fl23|) is valid precisely in orders O(e^) 
and O(e^), corrections of order O(e^) are nonlinear w.r.t. 6u- and (5f-shifts of the 
points A, B and C. 

Using ([23]) we calculate 5m(/i2) = ^'"(/oi) + ^w(/ii) - ^m(/o2) + ^3(^0, vq + e)e^ + 
o{e^) = 2fis{uo,vo)e^ + o{e^), ^^(/^a) = 6v^foi) - 6v{fl,) + 6v{fo2) + e{6v{n,) - 
Sv{foi))ui{uo,Vo + e) + ei5u{fo2)-Su{foi))iy2iuo,Vo + e) + iy3{uo,Vo + e)e^ + iy^{uo,Vo + 
e)e^ + o(e^) = e'^{dyiys{uo,Vo) + iyi{uo,vo)iys{uo,Vo)) + o(e^). As we see now, for 
points f^2 (respectively /21) only the shifts along the respective curvature lines is 
of order 0{e^) (only for surfaces of special type this shift may degenerate to 0), 
while the perpendicular shift is of order O(e^). For we now obtain a remarkable 
result: both its shifts are of order O(e^), compared to the error estimates of order 
0{e^) for Su{fl2) = 2€^{dufi3iuo,Vo) + ^2(^*0, ^'o)yU3(^io, fo)) + o{e^), 5v{f^2) = 
2e'^{dyU3{uo,vo) + i^i{uo,vo)u3{uo,vo)) + o{e^). 
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This observation suggests us to partition the complete discrete lattice //^ into 3 
sublattices: 

a) the even sublattice of points /2i,2j' 

b) the odd sublattice of points /2i+i.2j+i! 

c) the intermediate sublattice of points /|j_|_i_2i; f2i,2j+i- 

Infinitesimally, in a (A^e)-neighborhood of the initial point /qq (with <^ 1/^); 
the even sublattice has both Sud- and 5f£)-shifts of order O(e^), for the odd sublattice 
they have order O(e^), for the intermediate one only the shifts along the respective 
curvature lines are of order O(e^), the perpendicular shifts being of order O(e^) 
again. 

Accumulation of errors for these sublattices is also different: linear accumulation 
of 0(e^)-errors for the odd sublattice, quadratic for the even sublattice and mixed 
for the intermediate sublattice: linear accumulation of 0(e^)-shifts along the respec- 
tive curvature lines and quadratic in perpendicular direction: 5ti(/|j2j) — ^'j'^^i 
^<f2i,2j) - ^ ■ J ■ e^ ^M(/2Vi,2i+i) - {i + j) ■ e^ '5^^(/2Vi,2i+i) - {i+ j) ■ e^ 
^u{fli,2j+i) ^ ? ■ J ■ e^ 5v{f^i,2j+i) - + j) ■ e^- 

These estimates are easily obtained using ( j23l) (see below the proof of Theo- 
rem [T5ll , for example the first one follows from: 

<^^(/2Vl),20+l)) = ^<fki+^),2j) + <^«(/2i,2{i+l)) - <^^(/2i,2j) + ^^4) 

2e^{dufi3 + /i2/i3) - 2e/i2(5M(/|i,20+i)) - ^w(/2i,2i)) + o(e^). 

Note that the estimates of Lemma do not require the initial three points A, B, C 
to lie near the grid points fij of the original smooth net and the distances d{A, B), 
d{A, C) need not to be equal, both just have to be of order 0(e). 

Theorem 15 For a smooth surface without umbilic points parametrized by curva- 
ture lines f : Q —>■ M.^, Q = {{i,v)\u'^ + < 1} and sufficiently small e > 0, 
there exists a discrete circular net f with all its points fl^ on f{fl), such that 
dif.,Jt,) = 0{e^ 
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Proof. We will give the details for the global estimate Su{f2i2j) < C* ■ z ■ j ■ e^, the 
other are proved in the same way. 

First define a function S{k,n) = |5^(/2/fc 2(n+i)) ~ '^^(/2fc,2n.)l- From fl2^ we see 
that 

S{k + l,n)< Sik, n) + 2e|^2|5(/c, n) + 2e^\dufi3 + fi^l^^l + o(e^) 
and using Lemma [T2] we immediately obtain 

S{k,n) < Ce^k. 

Now for the functions Afc(n) = |5m(/|^, 2„)| taken separately for each column of the 
even sublattice one gets from (!24|) that 

Afc(n + 1) < Afc(n) + (1 + 2e|/i2|)^(fc, n) + 2e^K < Ak{n) + e^{2Ck + K) 
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with K = 2 max(„ t,)gf2 \dufi3 + Ai2/U3|- So fl22l) gives us the global estimate Ak{n) < 
C-k-n-e^n 

In [HI one can find a similar result with the same order of approximation but 
without the requirement f^j G f{^). 
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Appendix 

Lemma 16 Let an e-neighborhood f : ^ of a point M on a smooth surface 
be given, and for the principal curvatures one has Kl + > at M. Then the 
intersection Iabc of a plane ttabc lyi''ng at -distance from the tangent plane ttm 
to f{^le) at M is a curve lying in -neighborhood of a quadric — the corresponding 
Dupin indicatrix of the surface /(fi) in the central point M of the e-patch. 

Proof. Taking the appropriate Cartesian coordinate system with the origin M we 
can approximate the chosen e-patch by the osculating paraboloid 'J' = {z = {Kix^ + 
K2y'^)/2} with error terms of order o(e^). We obtain the quadric Qm = ^ H tiabc 
called the Dupin indicatrix of f{^) at M. As one can easily check, the angle between 
ttm and any of the tangent planes, taken at a point of f{fie), e^-close to Qm, is — e: 
the normals n = fu x fv to f{^e) at such points are e^-close to the normals ni of the 
osculating paraboloid at the corresponding points with the same (x, ?/)-coordinates. 
For the latter one has ni{x, y) = —Kixi + + k and for x ^ e and/or y — e, the 
angle between ni{x, y) and rii(0, 0) = A; is ~ e. 

Now using the fact that f{fle) is e^-close to the osculating paraboloid and stan- 
dard estimates for the values of implicit functions and their derivatives we conclude 
that Iabc and Qm are e^-close and the tangent directions at their e^-close points 
are also e^-close. □ 
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